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ABSTRACT 

The results of a three-dimensional model for disk-halo interaction are presented here. 
The model considers explicitly the input of energy and mass by isolated and correlated 
supernovae in the disk. Once disrupted by the explosions, the disk never returns to 
its initial state. Instead it approaches a state where a thin HI disk is formed in the 
Galactic plane overlayed by thick HI and HII gas disk with scale heights of 500 pc 
and of 1 to 1.5 kpc, respectively. The upper parts of the thick HII disk (the diffuse 
ionized medium) act as a disk-halo interface and its formation and stability are directly 
correlated to the supernova rate per unit area in the simulated disk. 
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; 1 INTRODUCTION 

This is the first of three papers where the dynamics of the 

■ Galactic fountain and its effects in the Galaxy are pre- 
' sented. The present paper describes three-dimensional hy- 
, drodynamical simulations of the disk gas powered by iso- 

■ lated and clustered supernovae and its interaction with the 
halo. 

In paper II (Avillez, 1999b) a discussion is conducted 
' on the dynamics of the Galactic fountain in the halo and the 
formation of HI clouds condensing from the fountain flow. 
In paper III a discussion of the present day modelling of the 
disk-halo interaction is carried out. This includes a detailed 
presentation and comparison between the existing models, 
their strengths and weaknesses and how their predictions 
compare with observations. 

The approach adopted here assumes the Galactic foun- 
tains are intimately related to the vertical structure of the 
thick gas disk and to the rate occurrence of supernovae per 
unit area in the disk. The hot gas released by the explosions 
breaks through the thick disk, escaping into the halo via the 
diffuse ionized medium, and forming large scale fountains. 

The term "galactic fountain" has been used to define a 
global circulation model where the disk gas rises up into the 
halo, cools and returns to the disk. The returning cold gas 
contitutes some of the observed HI clouds in the halo (see 
Wakker & van Woerden, 1997). The latest developments on 
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the study of the physics of the HI gas in the Galactic halo 
are reviewed in Wakker et al. (1999). 

Several mechanisms have been proposed for the disk- 
halo-disk circulation, taking into account the vertical struc- 
ture of the Galactic disk and the distribution of O and B 
stars in the disk. These comprise two major models: the 
classical fountain (CF) model (Shapiro & Field, 1976; Breg- 
man, 1980; Kahn, 1981) and the chimney model (Tomisaka 
& Ikeuchi, 1986; Norman & Ikeuchi, 1989). 

The CF models assumed the fountain gas being origi- 
nated in type II supernovae randomly distributed in the disk. 
As a consequence the gas gains enough energy to expand 
buoyantly into the halo where it cools and condenses into 
clouds that would rain over the disk. The maximum velocity 
these clouds would have is 70 km/s (Kahn, 1981), ruling out 
HVCs (clouds with descending velocities greater than 100 
km/s) as a product of a Galactic fountain. These models did 
not take into account the presence of a thick disk composed 
of a multiphase medium with different scale heights com- 
monly known as the Lockman (with a scale height ~ 500 
pc) (Lockman, 1984; Lockman et al, 1986) and the Reynolds 
(with a scale height /iz ~ 1 — 1.5 kpc) (Reynolds, 1987) lay- 
ers. The presence of such a thick gas disk constrains the 
dynamics of the hot intercloud medium on its way out of 
the disk. 

The detection of holes with sizes varying between a few 
hundred and a kilo parsec in the emission maps of nearby 
spiral galaxies (M31 and M33) by Brinks & Bajaja (1986) 
and Deul & Hartog (1990) and in the Milky Way by Heiles 
(1984) and Koo et al (1992) suggest that associations of O 
and B stars exploding in a correlated fashion could be re- 
sponsible for the formation of superbubbles that blow holes 



2 M. A. de Avillez 



in the disk. As a consequence, they release their hot inner 
parts into the lower halo through tunnels physically con- 
nected to the underlaying OB associations resembling chim- 
neys (Tomisaka & Ikeuchi (1986), Norman & Ikeuchi (1989), 
Mac Low et al, 1989). 

The energies involved in such a phenomenon are of the 
order of 10^^ erg (Heiles, 1991). As the hot gas rises up to a 
height of 10 kpc, it cools and returns to the Galactic disk, 
forming a fountain. The returning clouds would then acquire 
high velocities. Furthermore, chimneys reduce the volume 
filling factor of the hot gas in the disk reducing it from 0.7 
predicted by the three-phase model of McKee & Ostriker 
(1977) to some 0.2 as predicted by Norman & Ikeuchi (1989). 

Chimneys have been regarded as the most efficient way 
to expand the hot disk gas into the halo and as indirectly 
being reponsible for the formation of HVCs (Kuntz & Danly, 
1996; see reviews by Wakker & van Woerden 1997). The ma- 
jor drawback is the number of OB associations required to 
produce chimneys in such a way as to guarantee a continu- 
ous injection of hot gas into the halo. 

Rosen et al. (1993) and Rosen & Bregman (1995) de- 
veloped an integrated model where two co-spatial fluids rep- 
resenting the stars and gas in the interstellar medium were 
used. These models reproduced the presence of a multiphase 
media with cool, warm and hot intermixed phases having 
mean scale heights compatible with those observed in the 
Galaxy. However, the structure and properties of the ISM 
changed according to the overall energy injection rate. The 
best range of supernova rates that allow the reproduction 
of the ISM stratification are closed to the observed in the 
Galaxy. The resulting volume filling factors (50% for the hot 
medium, 25 — 30% for the warm gas and 20% for the cold 
gas (Rosen et al., 1999)) arc slightly different from those 
expected for the Galactic disk (see Spitzer, 1990; Ferriere, 
1995). 

Substantial progress has been made in modelling the 
disk-halo interaction and also the Galactic fountain since it 
was first suggested by Shapiro & Field in 1976. The models 
described above give a two-dimensional description of the 
dynamics of the disk and halo gas, and simulate the vertical 
distribution of IVCs and HVCs in the halo. Simplified con- 
ditions for the stars and ISM have been taken into account. 

The Galactic fountain models assumed that supernovae 
were randomly distributed in the disk and occurred with a 
rate of 3 per century, corresponding to a massflux of hot 
(T ~ 10^ K) material into the halo of the order of 10"^'' 
g cm s~^, whereas global models such as those of Rosen & 
Bregman (1995) showed that it is very difficult to reproduce, 
in the same simulation, both the structure of the ISM near 
the disk and, at the same time, the presence of HI gas with 
high and intermediate velocities in the halo without increas- 
ing the volume filling factor of the hot gas in the disk. 

Both the CF and the global models of the disk-halo in- 
teraction show that large scale outflows may be responsible 
for the formation of some of the observed IVCs and HVCs, 
providing the hot intercloud medium has a temperature of 
To = 10® K and the density p = 1.67 x IQ-^* g cm'^ 
This corresponds to mass influx of HVCs onto the disk of 
2 — 4 Mq yr^^ . A result in accordance with the estimates car- 
ried out by Wakker (1990). The chimney models neglected 
the presence of a dynamical thick disk which would provide 
constraints to the structure of the chimneys. The formation 



of the chimneys was carried out in a medium where no other 
phenomenon was present. 

Furthermore, the two-dimensional calculations impose 
natural limitations on the evolution of the gas and clouds. 
The absence of a third dimension constrains the motion to 
a vertical plane perpendicular to the Galactic disk. If three- 
dimensional evolution is considered, the overall structure of 
the flow may suffer modifications that lead to the generation 
of phenomena that may not be identified in two-dimensions. 

The principal objective of the research reported here has 
been to develop a three-dimensional model that can account 
for the collective effects of type lb, Ic and II supernovae 
on the structure of the interstellar medium in the galactic 
disk and halo. This model should, accordingly, account for 
the formation of the major features already observed in the 
Gala:x;y; these include the Lockman and Reynolds layers, 
large scale outflows, chimneys, and HI clouds, features that 
all contribute to the overall disk-halo-disk cycle known as 
the Galactic fountain. 

Section 2 deals with numerical modelling model where 
the model of the Galactic disk, rate of supernovae , boundary 
and initial conditions and numerical scheme are presented. 
Section 3 presents the evolution of the simulations, compris- 
ing the global evolution, the description of the simulations 
as they reach the steady state and the distribution of cold 
gas in the simulated disk. In section 4, a discussion on the 
simulations and the major predictions is carried out, and 
a comparison with existing observations is made. Finally, 
section 5 presents the conclusions and future prospects. 



2 NUMERICAL MODELLING 
2.1 Model of the Galaxy 

The study of the evolution of the Galactic disk rests in the 
realization that the Milky Way has a thin and a thick disk 
of gas in addition to a stellar disk. The thin gcis disk hcis 
a characteristic thickness comparable to that of the stellar 
disk of Population I stars. The thick gas disk is composed of 
warm neutral and ionized gases with different scale heights - 
500 pc (Lockman et al., 1986) and 950 pc (Reynolds, 1987), 
respectively. 

The stellar disk has a half thickness of 100 pc and the 

vortical mass distribution, p*, inferred from the star kine- 
matics, given by (Avillez et al, 1998), 

p* = p*,osec/i' [(27rG/3*p*)'/' z] (1) 

where z varies between —100 pc and 100 pc, p*,o = 3.0 x 
10^^'* g cm^'' is the mass density contributed by Popula- 
tion I stars near the Galactic plane (Allen, 1991) and the 
constant /9* = 1.9 x 10~^^ cm~^ s^. This mass distribution 
generates a local gravitational potential, of the form 

* = --^ In cosh [(27rGp*/3*)'/^ z] . (2) 

The stellar disk is populated by supernovae types lb, Ic 
and II, whose major fraction occurs along the spiral arms of 
the Galaxy close to or in HII regions (Porter & Filippenko, 
1987). Supernovae types lb and Ic have progenitors with 
masses M > 15 Mq , whereas Type II SNe originate from 
early B-type precursors with 9Mq < M < 15 Mq. 

The rates of occurrence of supernovae types Ib-|-Ic and 
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Volume 


Distribution 


(10-3 yr-l) 


(10-3 yr-l) 


(10-3 yr-l) 


Galaxy 


OB Assoc. 


1.2 


7.2 


8.4 




Isolated 


0.8 


4.8 


5.6 


Sim. Disk 


OB Assoc. 


0.002 


0.01 


0.012 




Isolated 


0.001 


0.007 


0.008 



Table 1. Rates of occurrence of supernovac types Ib+Ic and II in the Galactic disk having a volume Vq and in the simulated stellaj: 
disk with volume V . The rates account for isolated as well as clustered supernovae. 



II in the Galaxy are 2 x 10-3 yr-^ and 1.2 x 10 ^ yr i re- 
spectively (Cappellaro et at, 1997). Similar rates have been 
found by Evans et al. (1989) in a survey of 748 Shap ley- Ames 
galaxies. The total rate of these supernovae in the Galaxy 
is 1.4 X IQ-^ yr~^. Sixty percent of those supernovae occur 
in OB associations, whereas the remaining 40% are isolated 
events (Cowie et al., 1979). 



2.2 Basic Equations and Numerical Methods 

The evolution of the disk gas is described by the equations 
of conservation of mass, momentum and energy: 



0; 



a(pv) 
dt 

d{pe) 
dt 



-h V (pvv) = -Vp - pV$; 



-|- V (pev) = — pVv — n A; 



(3) 



(4) 



(5) 



where p, p, e, and v are the mass density, pressure and 
specific energy and velocity of the gas, respectively. The set 
of equations is complete with the equation of state of ideal 
gases. 

A is the functional approximation of the cooling func- 
tions of Dalgarno & McCray (1972), Raymond et al. (1976) 
and the isoclioric curves of Shapiro & Moore (1976), except 
in the range of temperatures of 10^ — 5 x 10^ K where the 
simple power law of the temperature (Kahn, 1976) 

A = 1.3 X io-i8r-°-^ erg cm3 (6) 

is applied. Between 5 x 10® and 5 x lO'^ the cooling function 



has been approximated by T 



(Dorfl, 1997). In order to 



avoid any cooling below zero temperature, A is zero below 
200 K. 

The equations of evolution are solved by means of a 
three-dimensional hydrodynamical scheme using the piece- 
wise parabolic method (Collela & Woodward, 1984) and 
the adapted mesh refinement algorithm of Berger & Collela 
(1989). Three levels of refinement are used. Each is created 
with a refinement factor of two. The maximum resolution 
obtained with this method is 1.25 pc. 



2.3 Computational Domain and boundary 

Conditions 

The simulations were carried out using a cartesian grid cen- 
tred on the Galactic plane with an area of 1 kpc^ and 
extending from -4 kpc to 4 kpc. Grid resolutions of 5 pc 



and 10 pc were used for —270 < z < 270 pc and for 
z < —250 pc and z > +250 pc, respectively. The cells lo- 
cated at -270 < z < -250 pc and 250 < z < 270 pc of 
the coarser and finer grids overlap to ensure continuity be- 
tween the two grids. The solution in the finer grid, at ranges 
-270 < 2 < -250 pc and 250 < z < 270 pc, results from 
a conservative interpolation from the data in the coarser 
cells to the fine grid cells as prescribed by Berger & Collela 
(1989). 

The boundary conditions along the vertical axis are pe- 
riodic; outflow boundary conditions are used in the upper 
and lower parts of the grid parallel to the Galactic plane. 



2.4 Initial Conditions 

The initial configuration of the system in the computational 
domain considers the vertical distribution of the disk and 
halo gases in accordance with the Dickey & Lockman (1990) 
profile, and reproduces the distribution of the thin and thick 
disk gas in the solar neighborhood. The gas is initially in 
hydrostatic equilibrium with the gravitational field, thus re- 
quiring an initial gas temperature as defined by the density 
distribution. 



2.5 Supernovae in the Simulated Disk 

Throughout the simulations the presence of the stellar disk 
is required with a thickness of 200 pc and an area of 1 kpc^, 
corresponding to a volume of V = 2 x 10* pc3. Isolated 
and clustered supernovae of types lb, Ic and II occur in the 
volume V at rates obtained from the values discussed in §2.1 
and given by 



Vg 



1.42 X 10"-' <T 



where Vg = 1.4 x IflU pc3 is the volume of the Galactic disk 
with a radius of 15 kpc and a thickness of 200 pc, and a is 
the rate of supernovae observed in the volume Vq. Table 1 
presents the rates of supernovae in the Galaxy as well as 
in the simulated disk. In the simulated stellar disk, isolated 
and clustered supernovae form at time intervals of 1.26 x 10^ 
yr and 8.4 x 10* yr respectively. 

The time interval between two sucessive superbubbles 
in the volume V varies between 1.9 x 10® yr and 5.23 x 10® yr 
(Fcrriore, 1995), but following Norman and Ikcuchi (1989), 
and adopting an average value of Ns at = 30 supernovae per 
superbubble, the time interval of formation of superbubbles 
in the volume V is 
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Figure 1. Grey-scale maps of the density distribution during the first 25 Myr for y = 500 pc at: (a) 1 Myr, (b) 5 Myr, (c) 10 Myr, (d) 
15 Myr, (e) 20 Myr and (d) 25 Myr after start of calculations. Note the collapsing of the gas during the first 15 Myr, forming a layer 
containing the major fraction of the initial mass. As the cold gas descends towards the plane it leaves behind a low density medium (light 
grey in images (b) and (c)). 



NsN „ _ ,„6 

— 2.5 X 10 yr, 



(8) 



the value used in this study. 

Each supernova is set up in the beginning of phase II, 
with a thermal energy content of (Kahn, 1975) 



Exhe 



and a kinetic energy content of 
Ekiu = 0.138poa, 



(9) 



(10) 



where a = 2E/po, po is the local density of the medium 
where the supernova occurs, and E is the energy of the ex- 
plosion. 

The generation of supernovae in the stellar disk is car- 
ried out in a semi-random way, where the coordinates of 
the supernova are determined from three new random num- 
bers converted to the dimensions of the stellar disk. After 
the determination of the random numbers, and therefore, 
the location of the new supernova in the grid, several tests 
are carried out in order to determine the type of supernova 
whether it occurs isolated in space or within an OB asso- 
ciation, or even to determine if the new location must be 
dismissed. 

The location of a new supernova in the grid is based on 
the following constraints: (a) vertical density distribution 
of population I stars given by equation (|l|), (b) type of new 



supernova, (c) rates of occurrence of supernovae types Ib-|-Ic 
and II and (d) rate of formation of superbubbles. 

When a supernova is due to appear in a specific location, 
a supernova generator checks for the occurrence of a previous 
supernova in that location. If this test is negative, then a new 
supernova is set up, providing its generation is compatible 
with the time interval between the previous supernova of 
the same type. If the test is positive, then the time delay 
between the old supernova and the new one is checked. Now 
the decision is based on the comparison of the this time delay 
and that of a new supernova within a superbubble. If the 
delay is smaller than the time delay within the superbubble, 
then no supernova is generated and the generator chooses a 
new position from a new set of random numbers, otherwise 
the supernova is set up. 

When the type of the new supernova is known, the algo- 
rithm determines the radius of the supernova in the begin- 
ning of phase II. This is acomplished by equating the amount 
of mass of the progenitor star released during its explosion 
and the mass of the ISM engulfed by the blast wave. The 
mass of the progenitor star is a ramdom value determined 
from the interval of masses characteristic to the progenitors 
of that type of supernova, e.g. if the next supernova to occur 
is a type lb, then the amount of mass that is released during 
the explosion is found within the interval [15, 30] Mq. 
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3 EVOLUTION OF THE SIMULATIONS 



The simulations were carried out for a period of 1 Gyr with 
supernovae being generated at time 0. The simulations start 
from a state of hydrostatic equilibrium breaking up during 
the first stages of the simulation. The system evolves into a 
statistical steady equilibrium where the overall structure of 
the ISM seems similar on the global scale. 



3.1 Global Evolution 

During the first 50 Myr of simulated time there is an im- 
balance between cooling and heating of the gas in the disk. 
Initially there is an excess of radiative cooling over heat- 
ing because of the small number of supernovae during this 
period. 

The gas originally located in the lower halo starts cool- 
ing and moves ballistically towards the midplane, colliding 
there with gas falling from the opposite side of the plane 
(Figure 0). After 15 Myr of evolution, the major fraction of 
the initial mass is confined to a slab having a characteristic 
thickness of 100-150 pc (Figure ^ (d)). A thin wiggly disk of 
cold gas forms at the midplane and has a thickness of few 
tens of parsec. As the supernovae occur, they warm up the 
gas in the slab, gaining enough energy to overcome the grav- 
itational pull of the disk and, therefore, expanding upwards 
(Figure |l (e) and (f)) redistributing matter and energy in 
the computational domain. 

After the over-expansion, a dynamic equilibrium is set 
up between upward and downward flowing gas (Figure |^. 
During the first 10-20 Myr, the descending gas has a max- 
imum value of 3.4 x 10"^^ g cm"^ s"^ (that is 4.8 x 10"^ 

The ascending gas massfiux peaks at 
2 7i/r 1. — — 2 



2.8 X 10"^^ g cm"^ s"^ (3.9 x 10"^ Mq kpc"^ yr"^). The 
expansion of the disk gas diminishes during the next 20 Myr, 
with the consequent decrease in the total massfiux. Between 
60 and 75 Myr, a balance between descending and ascending 
flows sets in. Both the ascending and descending gases have 
massfiuxes of approximately 3x 10~^° g cm~^ s~^ (4.2 x 10"'^ 
Mq kpc~^ yr~^), which corresponds to a total infiow rate 
of 2.97 Mq yr~^ in the Galaxy on one side of the Galactic 
plane. 

Such a dynamic equilibrium is observed during the rest 
of the simulations, indicating that a quasi-equilibrium be- 
tween the descending and ascending gases dominates the 
evolution of the disk-halo gas. However, there are some pe- 
riods when one of the fiows dominates over the other. This is 
related to periodic imbalances between cooling and heating 
leading to compressions and expansions of the disk gas. 

As the balance between descending and ascending flows 
is reached, the disk evolves into a statistical steady equilib- 
rium where, in general, the disk gas is composed of a mul- 
tiphase medium extending into the halo, forming a frothy 
region laced by sheets and crossed by HI clouds and well 
collimated chimneys as can be seen in Figure ^ The figure 
presents the vertical distribution of the density up to 1 kpc 
on either side of the plane. 

The disk gas shows a stratified distribution where 
neutral and ionized phases co-exist having different scale- 
heights. The neutral phase is composed by cold (T < 10^ 
K) and warm (10^ < T < 10* K) gas whereas the ionized 
phase is composed by a warmer gas with 10* < r < 10^ 
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Figure 2. Total massfiux cj> = pvz of the ascending (vz > 0, 
shown as diamonds) and descending {v^ < 0, shown as squares) 
flows measured at z = 140 pc and averaged over the entire area of 
the Galactic disk during the first 500 Myr of evolution. A balance 
between the ascending and descending flows sets in after 75 Myr. 



K (hereafter identified as warm ionized gas) and hot gas ( 
T > 10^ K). 

The neutral component is mainly distributed in the 
layer \z\ < 500 pc, forming a thick HI disk and having an 
average temperature of 10* K. Cold gas is present, predom- 
inantly in a thin but irregular layer (thin HI disk) around 
the Galactic plane, having a characteristic thickness of tens 
of parsec (Figures]^ and^, whereas the hot gas is mainly 
distributed above 1.5 kpc. 

Gas with temperatures varying between 10* and 10^ K 
fill up to ~ 80% of the disk volume, whereas gas with tem- 
peratures smaller than 10* K, fills on average (Figure 
50% of the disk volume. The T < 10* gas is mainly located 
in a layer with half thickness of 500 pc {\z\ < 500 pc) (see 
also Figure ^), although some of this gas extends to 1.5 kpc, 
but with a very small volume filling factor. The warm ion- 
ized gas dominates the region between 500 pc and 1.5 kpc, 
although its presence is detected up to 2.5 kpc. Gas with 
temperatures larger than 10^ K fill on average 15 — 23% 
of the disk volume below 1 kpc. With the increase of z, 
the volume occupied by the hot gas increases, becoming the 
dominant medium at 2 > 1.5 kpc (Figure |^. The density 
distribution of the disk gas decreases with z. At \z\ < 200 
pc it can be approximated by a gaussian followed by an ex- 
ponential distribution (Figure ^). The gaussian is overshot 
by the exponential distribution with a smooth decrease fol- 
lowed by a steep increase in the density profile (Figure ^ 
on either side of the plane. This suggests the presence of an 
interface between the thin disk and the thick HI disk. Such 
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X (pc) X (pc) 

Figure 3. Color maps showing the structure of the Galactic disk measured at two locations in the disk: y = 480 ((a) and (b)) and 
y = 800 pc ((c) and (d)). The times at which these images refer are: (a) 793 Myr, (b) 794 Myr, (c) 982 Myr and (d) 983 Myr. The color 
scale varies from dark blue (low density gas) to red (high density gas). Neutral gas is represented by a range of colors varying between 
yellow and red, whereas ionized gas is shown in blue scale. The images show that globally the disk has similar structures, although 
locally there are variations. Note the global expansion of warm (light blue) and hot (dark blue) ionized gas carving the thick disk and 
the constant presence of the wiggly cold disk (red) at 2 = pc. Local outbursts such as mini-chimneys ((a) and (b), located at x = 500, 
2 = —100 pc) and chimneys ((c) and (d)) are also observed. 
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an interface is fed by the hot gas flowing upwards from the 
thin disk in large scale outflows {I ~ 200 pc or more) inter- 
acting with the cooler thick disk gas overlaying it. Such a 
configuration becomes Rayleigh- Taylor unstable, and there- 
fore the upward flow acquires a finger-like structure during 
its expansion through the thick disk. 

This density distribution gives the appearance of a 
stratified distribution which is compatible to the volume fill- 
ing factors of the difi'erent phases of the ISM at different z 
(see Figure ^). The ionized component of the disk gas ex- 
tends upwards up to z ~ ±1.5 kpc. Its density decreases 
smoothly up to 1.4-1.5 kpc, where it has a steep decrease. 
This region acts as an interface between the thick disk and 
the halo gas. For z > 1.5 kpc the gas has lower densities and 
temperatures T > 10^ K. 



3.2 Statistical Steady Evolution 

The disk gas shows a steady equilibrium where its struc- 
ture and dynamics are similar throughout the rest of the 
simulations. That is, the disk has the same structure at dif- 
ferent times and locations, although it experiences changes 
due to local events such as supernovae, superbubbles, chim- 
neys, collisions with HI clouds, shock waves, etc. In order 
to stress this point, images showing the time evolution of 
the disk gas at different locations are presented in Figures 
^|ll|. The figures were taken in directions perpendicular and 
parallel to the midplane at 500 and 800 Myr of evolution 
time. The a;z— slices (Figures were taken at 30, 220, 410, 
600, 790 and 980 pc from the edge of the computational do- 
main, whereas the xy— slices (Figures p^pT[ ) were taken at 
z = —250, and 250 pc. The density distributions of the 
disk gas shown in Figures |^, ^ and ^ are complemented 
with the corresponding temperature maps in Figures |^, ^ 
and|l^. 

The figures show a large number of features comprising, 
among others, HI clouds and sheets, bubbles and superbub- 
bles, interacting supernovae forming networks of hot gas, 
large regions of hot gas crossing the arj/— slices at z — —250 
and 250 pc, etc. These features can be summarised as fol- 
lows: 

• A wiggly thin and a cold disk, having a characteristic 
thickness varying between a few parsec and tens of parsec 
(Figures There are regions where the thin disk may 
have disappeared as result of a large number of supernovae in 
the neighborhood as shown in Figure ^ (a) at 200 < x < 400 
pc. The region is crossed by a stream of hot gas as can be 
seen in the corresponding temperature map. 

• Vertical structures connected to the underlying thin 
disk and appearing to wiggle around like worms. These 
structures are the debris of broken shells and supershells, 
that occurr displaced from the midplane. 

• Sheets of cold gas distributed in the disk that result 
from the breakup of larger structures (Figures ^ - ^l|) . The 
sheets have thicknesses of a few pc and widths of several 
tens of parsec, even hundreds of parsec. The sheets observed 
at z = ±250 pc seem to connect to each other, forming a 
"network" of cold gas. In the plane, the cold gas forms large 
regions of very high density and low temperature (Figure hd 
(b) the darkest blue regions). 

• Small cloudlets in the vicinity of worms and other sheet- 
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Figure 6. Gray scale maps of the gas density measure at 500 
Myr of disk evolution for; (a) y=30, (b) 220, (c) 410, (d) 600, (e) 
790 and (f) 980 pc from the edge of the computational domain. 
A thin HI disk with a variable thickness and a wiggly structure 
is present in all the images. Its number density varies between 1 
and 50 cm"''. The largest density values are found in the mid- 
plane stretching over large regions with sizes of 50 — 150 pc. Cold 
gas is mainly distributed in the thin disk and in sheet like struc- 
tures found isolated or connected to the underlying thin HI disk, 
resembling worms. 
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Figure 4. Density profile of tlie disk and lialo gas measured along two line of sights passing through the plane at: (a) x = 300, y = 250 
pc and (b) x = 750, y = 750 pc. Note the presence of a thin disk with a thickness of 100 pc centred at z = 0. The density profile behaves 
like a gaussian distribution up to 200 pc on either side of the plane followed by an exponential up to 1 kpc. Above this height the density 
decreases steeply showing a separation between the warm ionized medium and the halo gas at 1.5 kpc. 
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Figure 5. Best fits for distribution in z of the volume filling factors of cold (solid line located at the lower left edge of the plots), warm 
neutral (dotted line), warm ionized (dot-dashed line) and hot (long dashed line) gas. The fits were calculated from a series of 50 profiles, 
measured with a time interval of 1 Myr, and taken at: (a) 450 — 500 Myr and (b) 950 — 1000 Myr. The variation of the volume filling 
factors with z are compatible with a stratified distribution of the disk and halo gases. The neutral gas dominates for z < 500 pc, the 
warm ionized gas is mainly found in the layer located between 500 pc and 1.5 kpc and the hot gas dominates for z > 1.5 kpc. The height 
at which the hot gas becomes dominant is identified as the disk-halo interface. 



like structures. The smallest cloudlets detected in the images 
have 1.25 pc size - this corresponds to the maximum resolu- 
tion used in the disk. The cloudlets result from the breakup 
of the worms and other clouds which have been swept up by 
blast waves, involved in cloud-cloud collisions, etc. 

• Supernovae and superbubbles at different phases of 
their evolution surrounded by shells with variable thick- 
nesses. Supernovae and superbubbles located above the thin 
gas disk show an elongated appearance as result of the den- 
sity gradient along the direction (Figures ^j^. 

• Networks of hot gas connecting supernovae. These can 
be observed in both the xz— slices and a;i/— slices of the tem- 
perature distribution of the disk gas. The network of hot 
gas is not confined to the plane {z = 0), but spreads in the 
volume of the disk between z — —250 and 250 pc. 

• Reservoirs of hot gas (large regions of hot gas) moving 
away from the plane (Figures^, ^ and pj] ). 



In general the disk gas is pervaded by isolated supernovae 
and superbubbles (large supernova structures generated by 
sucessive supernovae). Supernovae change the local struc- 
ture of the interstellar medium, but are unable to change 
the global structure. Individual and clustered supernovae 
dominate the local environment, depending on their spatial 
location. Isolated supernovae change the structure of the in- 
ner parts of the disk whereas supernovae in OB associations 
dominate the upper regions of the disk. 

As the blast waves from supernovae expand in the disk, 
they sweep up clouds and sheets, triggering their disruption 
into smaller structures. At 2 = ±250 pc the presence of 
sheets results from: (a) cold gas descending from above and 
acquiring sheet-like forms as a result of Ray leigh- Taylor in- 
stabilities that occurr in larger clouds and (b) the breaking 
up of shells and supershells that expand upwards and are 
displaced from the Galactic plane. 
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Figure 7. Temperature maps of the same regions shown in the 
last figure. Note the presence of the thin HI disk with tempera- 
tures varying between 10^ and lO'^'* K. The coldest regions are 
those with larger densities (compare with previous figure). The 
major fraction of the disk is filled with gas having temperatures 
of some 10* K. Hot gas is found mainly in regions where super- 
novae went off and constitutes a network of hot tunnels that may 
cross the midplane. 



^1 



Figure 8. Gray scale maps of the gas density measure at 800 
Myr of disk evolution for; (a) y=30, (b) 220, (c) 410, (d) 600, (e) 
790 and (f) 980 pc from the edge of the computational domain. 
These images are similar to those shown in Figure ^ although 
locally there are changes due to different events occurring in the 
disk. 
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Figure 9. Temperature maps of the same regions shown in the 
last figure. Note the presence of streams of hot gas crossing the 
midplane, but this is a locallized event resulting from the large 
number of explosions that occurred in this region in a small period 
of time. 




zoo 400 600 BOO 1000 

X (pc) 

,^ - 500. Myr (h) 



i 



200 400 600 BOO 1000 




200 400 



Figure 10. Maps of the density distribution in the disk at 500 
Myr and for: (a) z = -250, (b) 2 = and (c) z = 250 pc. 
The HI gas observed ai z = ±250 pc is mainly concentrated in 
sheet-like structures, connected to others, some of the sheet are 
the walls of shells and supershells. Some of the cold gas present 
at 2 = is distributed in large regions whose centres have the 
largest densities - the darkest blue regions in image b). 



Young supernovae expand into a highly unstructured 
medium interacting with other supernovae or cold gas in the 
form of sheets or clouds (Figures ^-^ij) . During their expan- 
sion, the remnants are preceeded by shock waves engulfing 
the shells of other supernovae and triggering their disrup- 
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Figure 11. Color maps of the disk temperature at 500 Myr and 
for: (a) z = -250, (h) z = and z = 250 pc. 



tion. If the average shock speed in the shell is t;^ ~ 1 x 10* 
km s~^, then the transmitted shock takes a time 



AR 



(11) 



to cross the old shell, AR being the thickness of the shell 
normally 5 to 10 pc thick. Thus, the blast wave crosses the 
shell in tcross < 10* yr, generating a large pressure gradient 
which is enough to disrupt the old shell. As a consequence, 
the hot gas flows freely between the remnants forming a 
network of channels of low density gas. On average, after 



5 X 10^ years, a type II supernova releases 295 Mq of hot 
gas into the surrounding medium. Therefore, when several 
isolated supernovae merge together, they contribute to the 
formation of "reservoirs" of hot gas with enough energy to 
overcome the gravitational pull of the stellar disk expanding 
buoyantly upwards (Figures ^, ^ - ^) ■ 

In its ascending motion, the hot gas interacts with the 
denser gas distributed in the thick disk. Such a configuration 
is Rayleigh- Taylor unstable and as a consequence, the hot 
gas acquires a finger-like structure with a mushroom cap. As 
it moves upwards it appears to carve the cooler layers that 
compose the thick disk (Figure ^ . 

Bubbles and superbubbles occurring in the outer parts 
of the thin disk, at some ~ 100 pc above and below the 
Galactic plane, expand vertically due to the local stratifi- 
cation of the ISM becoming elliptical (Figures ^ and 
The vertical elongation of the shell surrounding the bub- 
bles/superbubbles depends not only on the local medium, 
but also on the amount of energy released by the super- 
nova explosions involved in the production of the bub- 
ble /superbubble. 

During their vertical expansion the shells, eventually 
become Rayleigh- Taylor unstable and blow their caps, re- 
leasing their hot inner parts into the ISM through narrow 
channels. The vertical expansion of the hot inner parts of the 
bubbles is halted by pressure effects due to the thick disk 
gas, while superbubbles break through it, blowing holes in 
the thick disk gas. The presence of such structures is shown 
in Figures | and |l§|l3[ 

These structures may be classified as mini-chimneys and 
chimneys depending on their vertical elongation and on the 
presence of actual tunnels breaking through the disk gas as 
can be seen in Figure ^ (c) and (d). The tunnel observed in 
the figure crosses the thick gas disk and has a characteristic 
width of approximately 120 pc. The tunnel maintains the 
same appearance providing enough supernovae occur in the 
underlying association. 

The expansion of hot gas in well coUimated structures 
may not have enough energy to break through the gas disk 
(Figure ^ (a) and (b)). Such a structure may be classified 
as a mini-chimney. The walls of the tunnels are formed by 
the remains of the broken shells whose appearance resem- 
bles worms crawling out of the disk. The average width of 
mini-chimneys varies between 50 and 60 pc whereas that of 
chimneys ranges from 100 pc to a few hundred parsec, de- 
pending on the number of supernova explosions inside the 
superbubble. Figure shows the sequential evolution of the 
disk gas in a region where a chimney is formed, while Figure 
presents the HI positional maps of the same region. The 
chimney, formed from a superbubble with a vertical elonga- 
tion of 200 pc and a width of 100 pc, is located in the region 
^ 350 < X < 450 pc in the southern hemisphere (hereafter 
called region A). The base of the superbubble is located in 
the midplane and connects with two merging bubbles lo- 
cated in the northern hemisphere at 300 < a; < 480 pc and 
40 < z < 150 pc. After a million years, these two bubbles 
merge, creating a large cavity extending from 350 < x < 480 
pc and 50 < z < 190 pc (Figures |l2| and |l^ (b)) and here- 
after named region B. 

At time 203 Myr, a supernova with radius of some 50 
pc (the supernova is well inside phase III of its evolution) is 
observed to occur at x — 400 pc, z = —100 pc, inside region 
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Figure 12. Evolution of a cliimney formed in an OB association located at {x = 400, y = 30, 50 < 2 < 150) pc. The images were taken 
at: (a) 203 Myr, (b) 204 Myr and (c) 205 Myr of simulated time. This figure is particularly rich in details regarding to the time evolution 
of the superbubble. After the superbubble has been formed sucessive explosions displaced the upper parts of its shell as can be seen in 
figure (a) located at z =~ —200 and —350 pc. The shell caps have been separated from the shell and conserve their forms. Note the 
increase in thickness of the left wall of the bubble as sucessive explosions occur. A large plume of hot gas rises near the superbubble at 
< X < 200. This structure expands as result of the high energy content in the reservoir of hot gas lying underneath. As the gas flows 
upwards it acquires a finger-like structure, whereas the chimney conserves its collimated tube-like form. 






Figure 13. HI emission maps of the previous figure. The figure presents twenty equally spaced contours with column densities varying 
between 10^* and lO'^" cm"'^ of the disk gas. One of the walls of the chimney is destroyed by the expansion of supernova located in its 
vicinity. The remaining wall has a thickness of several parsec and an elongation of 220 pc. At {x = 150, z = 200) pc a shell is moving 
upwards (image (a)) detaching from other cold gas located at {x = 220, z = 150) pc (image (b)) and later it break up into small clouds 
with sizes of several parsec (image (c) at 100 < x < 200, 150 < z < 240 pc). 



A. similar explosions occurred in the superbubble leading to 
its expansion and break up of the top of its shell (Figure |l^ 
(b) . This is corroborated by the HI emission map in Figure 
|l3| (b). A narrow tunnel of low density gas is surrounded 
by thick walls (thicknesses of a few parsec). The top of the 
chimney has a width of some 80-90 pc. 

The hot inner parts of the superbubble are released up- 
wards through sucessive injections due to different super- 
novae occurrences inside the cavity. As the gas expands out 
of the superbubble, it pushes upwards from z — —300 to 
z = —450 pc, debris of old shells in a period of 200 Myr. 
Continuing explosions inside the superbubble provoke the 
escaping of hot gas not only in the southern hemisphere, 
but also inject gas into region B (Figure |l^ (c)). As a con- 
sequence the upper parts of the shell surrounding region B 
also breaks up and an injection of hot gas into the thick disk 
through a mini-chimney is induced. 

The gas released by the mini-chimney is unable to break 



through the thick disk and acquires a finger-like configura- 
tion until it cools down to a recombination temperature and 
merges together with the surrounding medium (Figure |l^ . 
The figure shows the evolution of the same region in space 
as Figures |l| and but at later times: 209-211 Myr. After 
a period of a few million years, the width of region A varies 
between 200 pc near the plane to 110 pc on the upper parts 
of the sheU (Figure |l| (c)). 

3.3 Cold Gas in the Disk 

HI emission maps showing the vertical distribution of the 
cold gas in the disk and the corresponding temperature 
maps, taken at along the plane crossing the midplane at 
y — 220 pc, are presented Figures |l^^ The figures present 
a time evolution of the cold gas in the disk during a period of 
four million years and show 20 equally spaced contours rep- 
resenting values of column density between \og{NHi) = 18 
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Figure 14. Time evolution of the temperature of tlie disk gas measured at j/ = 220 between 654 and 657 Myr. The brightest regions 
are cold gas whereas the hottest regions are dark grey. The thin HI disk has a temperature of some 10^'* = 630 K, whereas the hottest 
regions have temperatures of ~ 10^'^ — 10^ K. The thin HI disk has a temperature of some lO'^'* = 630 K and a thickness ranging from 
tens of pc to a hundred parsec or more in some locations. The hot region located at (280 < x < 350, 70 < z < 150) pc (image (a)) 
expands acquiring a mushroom-like structure as result of the interaction with the cooler medium above it. The mushroom is already 
formed at 655 Myr (image (b)) and acquires an increased elongation as it evolves (images (c) and (d)). A stream of hot gas is observed 
at the base of the mushroom which is connected to the underlying hot source. 
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Figure 15. Emission maps of the disk gas shown in Figure The figure presents ten equally spaced contours representing column 
densities log{NHi) = 18 and log{NHi) = 20. The images show the presence of small HI clouds with a filamentary structure and connected 
to the underlying disk. Supernova cavities are surrounded by shells whose column density decreases with z. 



and log(A''£f/) = 20, where Nhi is in cm~^. The concen- 
tration of contours in a specific direction indicates a large 
column density of the gas there. 

In general the figures show HI emission concentrated in 
a wiggly cold structure (thin HI disk) distributed along the 
Galactic plane, having a thickness of a few tens of parsec, 
and in some cases being as thin as 10-20 pc. The presence 
of a very thin HI disk (with Az ~ 10 — 20 pc) is an occa- 
sional phenomenon and results from local variations in the 
structure of the thin disk as result of, for example, corre- 



lated supernovae located near the plane or collisions of high 
velocity clouds with the thin disk. The occurrence of the su- 
pernovae leads to the formation of tunnels of hot gas crossing 
the Galactic plane, and therefore, to the absence of the thin 
disk observed in some of the a;2:-maps (see Figure 9). Colli- 
sions of high velocity clouds, may compress and drag with 
the cloud parts of the thin disk. The largest column density 
observed in the images is caused by the low scale height of 
the cold gas. 

The figures show the presence of low column density re- 
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Figure 16. HI emission of the cold gas at z=0 pc and taken at: (a) 251 Myr, (b) 252 Myr, (c) 253 Myr and (d) 254 Myr of simulated 
time. The images show the time evolution of the cold gas in the midplane. The images show a sheet of cold gas with holes and channels 
of low density material forming a network of hot gas. Small clouds with sizes of few parsec result from the break up of larger structures. 
The presence of holes is a transient phenomenon that depends on the location of new supernovae as well as from the time of cooling of 
the disk gas. 



gions (absence of contours) persisting for long periods. These 
regions result from supernovae that occurred some time in 
the past, and that have blown holes in the Galactic disk. As 
the expanding gas cools, a thin shell of colder gas is formed 
behind the blast-wave generated by the supernova. 

The vertical structure of the HI holes is observed in 
most of the figures. Some are partially surrounded by thin 
sheets of cold gas, which result from old shells. The cone- like 
structures of low column density surrounded by cold sheets 
are found located on each side of the Galactic plane. Because 
of the high column density in the plane, the cones slightly 
change its local structure, but are unable to break the thin 
disk. 

Shells displaced from the midplane have a non-spherical 
appearance, resulting from the local medium not being 
smooth (on the small scale - 5 pc) and partly from the steep 
variation of the density along the z— direction - the shells 
acquire a ovoid form (Figures |l3[ |l^ - |l5| ). As the bubbles 
expand upwards the shells break, releasing the hot inner 
parts of the remnant into the surrounding medium. There- 



fore, the remaining parts of the shells become the walls of the 
cone-like structures. These walls have a wiggly appearance 
resembling "worms" crawling into the halo. 

Face-on maps of HI emission from the midplane, taken 
in a six million years period (Figure |l^), show high column 
density regions spread over the plane surrounding holes and 
channels of low column density gas. The holes and tunnels 
form a large scale network of low density gas spread over the 
Galactic disk. The tunnels have widths of some 100 pc and 
lengths in the range of 200-400 pc. This network forms as a 
result of young supernovae occurring near to old remnants. 
High velocity shock waves sweep up the old remnants, dis- 
rupting their shells and leading to a movement of the hot 
gas inside the young remnants into the low density medium 
surrounding it. The blast wave may suffer several reflections 
inside the tunnel, leading to the formation of small clouds as 
result of the disruption of the shells and condensation into 
clouds of the gas swept several times by the blast waves. 
This gas is compressed during each swept up to four times 
the local density (assuming that these shocks are strong). 
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Figure 17. Column density cm ^) of cold gas in the disk as a function of the velocity (km s ^) along one column (taken 

at —100 < z < 100 pc). The extreme values of the FWHM are 6.1 km/s and 10.5 km/s with the average values varying between 9.3 and 
9.8 km/s. 



The increase in the density leads to an increase in the cool- 
ing of the material and therefore to its condensation into 
clouds. 

Most of the gas observed in the HI emission maps 
forms clouds with sheet like appearance and connected to 
each other. Clouds have sizes/thicknesses varying between 
few parsec to hundreds of parsec. The larger clouds suffer 
Rayleigh- Taylor instabilities and break up into sheet-like 
structures that fall towards the Galactic plane. However, 
isolated structures with large sizes, identified near the plane 
are moving towards it colliding with it. The collision leads 
to the compression of the thin disk where higher density gas 
accumulates forming possible giant molecular clouds where 
star formation would occur. 

Column density profiles of the HI gas along a column 
of 400 pc (with —200 < z < 200 pc) versus velocity are pre- 
sented in Figure The figure presents profiles convolved 
with a Gaussian that has a velocity dispersion that corre- 
sponds to the local gas temperature, being each cell of the 
column normalized to the local gas density. The profiles are 
symmetric with respect to the central peak and have FWHM 
values varying between 6.1 and 10.5 km s~^, being these the 
extreme cases. The average values of FWHM are ~ 9.5 km/s. 

The major fraction of cold gas in the disk has positive 
and negative velocities varying between —20 and 20 km/s. 
However, there are clouds with larger velocities indicating a 
fast approach to the midplane, where they collide with cold 
gas falling from the opposite hemisphere, or an escape from 
the thin disk (positive velocities). 



4 DISCUSSION 

The results of large-scale modelling of a section of a Galactic 
arm, located at 8.5 kpc from the Galactic centre are pre- 
sented in this paper. It is a three-dimensional model that 
assumes the ISM as an integrated system where the collec- 
tive effects of supernovae (isolated and clustered) are taken 
into account. The supernovae occur in a medium strongly 
structured by previous supernovae, ionization fronts and ra- 
diative cooling. 

The model considers the interstellar gas distributed in 
a smooth thin disk with the vertical distribution of the cool 
and warm neutral HI gas given by Dickey & Lockman (1990) 
and warm ionized gas given by Reynolds (1987). It expli- 
citly includes both correlated and random supernovae in a 
manner compatible with observations. Sixty percent of the 
supernovae occur within associations, whereas the remaining 
occur isolated. The supernovae belong to types lb, Ic and II, 
and their masses are taken into account, as mass loading is 
allowed. 

Once disrupted by the explosions, the disk never returns 
to its initial state. Instead, it approaches a state where a thin 
HI disk is formed in the Galactic plane, overlayed by a thick 
disk of warm gas. 

The thick disk is generated as long as enough super- 
novae occur in the disk, regardless of the initial density 
profile adopted for the disk gas (§3.1), and, therefore, its 
formation and stability are directly correlated to the super- 
nova rate per unit area in the simulated disk. The density 
distribution of the disk gas in z can be approximated by a 
linear combination of two gaussians and an exponential with 
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Figure 18. Best fits (solid line) for n{HI) distribution in z (dashed line) along four lines of sight crossing the disk plane at: (a) 
{x = 30, 3/ = 30) pc, (b) (x = 300, y = 250) pc, (c) (x = 600, y = 500) pc and (d) (x = 800, y = 800) pc. Note the presence of humps in 
the n{z) distribution on each side of the thin disk. These humps prevent the hot gas coming from the thin disk to expand freely into the 
thick disk. 



Distribution A B C D E F Hz 

(cm~^) (pc) (cm~'^) (pc) (cm"'') (pc) (pc) 



Fig. 
Fig. 
Fig. 
Fig. 



18 (a) 
U (b) 

y (c) 
li (d) 



0.395 
0.90 
0.81 
0.70 
0.76 



90 
30 
50 
50 
54 



0.107 225 

0.110 200 

0.101 200 

0.105 205 

0.105 208 



0.064 
0.07 
0.04 
0.08 
0.09 



403 177.6 

405 158.6 

405 191.7 

405 182.4 

405 203.6 



Table 2. Values of the coefficients A, B, C, D, E and F in equation (|l^) and effective scale height Hz for the density profiles of Dickey 
& Lockman (1990) and of Figure Fsl (a)- (d). 



the general form, based on the model of Dickey & Lockman 
(1990) 



n{z) 



0.566 



Ae 



(12) 



where n is the number density of the ISM phase, Uo is 
the number density at the Galactic plane. A, C and E axe 
paremeters written in units cm~^ and B, D and F are writ- 
ten in pc. The density distribution yields an effective scale- 
height 



Noting that 



n{z)dz. 



-z2 0r 
e dZ — - — , 
2 ' 



(13) 



(14) 



and using (^3|) the effective scale-height of the density dis- 
tribution in equation ( p^ is given by 



Hz = 



(AB + CD)^/- 



EF 



(15) 



0.566 

For the density distribution of the disk gas in the Galaxy, 
Dickey & Lockman (1990, hereafter denoted by the DL dis- 
tribution) inferred the parameters shown in the first row 
of Table 3. The table presents, in rows 2-5, the parameters 
corresponding to density distributions measured along four 
lines of sight crossing the disk at: (a) x = 30, y — 30 pc, 
(b) X = 300, y = 250 pc, (c) x = 600, y = 500 pc and (d) 
X = 800, y = 800 pc. These distributions are the best fits to 
the underlying curves of the density distribution as shown 
in Figure 

As can be seen in Table 3, the coefficients C,D,E and 
F are similar between all the fits and the DL distribution, 
whereas the A and B parameters for the fits and for the DL 
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Component 
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H, 




(cm^ 




(K) 


(kpc) 


Thin HI disk 


0.6 - 


1 


< 8 X 10^ 


0.09 - 0.12 


Thick HI disk 


0.1 




10" 


0.5 


Thick HII disk 


0.01 - 


0.02 


10** - 10-^ 


1 - 1.5 



Table 3. The table presents the physical properties of the thin 
and thick disk gas in the simulated volume, n, T and Hz denote 
the volume density, temperature and height of the disk gas. 



distribution have variations as much as 15%. The diflFerence 
in the extreme case indicates that, in the simulations, most 
of the cold gas in the disk is concentrated near the plane 
and thus has a smaller effective scale-height. 

Despite the differences between the DL distribution and 
those obtained in the simulations, the thick gas disk has a 
distribution compatible with the presence of two gas com- 
ponents having different scale heights: a neutral layer with 
a scale height of 500 pc (warm HI disk) and an ionized com- 
ponent extending to a height of 1 to 1.5 kpc above the thin 
HI disk. 

Similar distribution of the thick gas disk and the HI and 
HII thick disks has been observed in the Galaxy (Lockman, 
1984; Lockman et al., 1986; Reynolds 1987) and in external 
galaxies such as NGC 891 (Dettmar, 1992). Studies of the z- 
distribution of HI in edge-on galaxies (NGC 4565 and NGC 
891) carried out by Rupen (1991) have shown the presence 
of an HI layer around the galactic midplane. This layer has a 
thickness that varies from 120 pc to 300 pc in both galaxies. 

The average physical properties of the different phases 
in the simulated disk in the steady state evolution are sum- 
marized in Table 3. These values are in agreement with those 
derived by Karberla & Kerp (1999) for the best-fit parame- 
ters required for a statistical steady equilibrium of the halo 
and disk gas in the Galaxy. 

4.1 Volume Filling Factors of ISM Phases 

The volume filling factors of the different ISM phases have 
been presented in §3.1 and the corresponding best-fit profiles 
are presented in Figure ^. The distribution of the hot gas 
varies from a occupation volume of some 23% in the inner 
disk (below 500 pc) to 100% above 2.5 kpc. 

The thin disk is populated by a cold medium with a 
volume-filling factor smaller than 5%. The warm neutral gas 
dominates below 500 pc occupying a volume of some 50%, 
whereas the warm ionized medium dominates between 500 
and ~ 1500 pc. At some 1500 pc, both the warm ionized and 
hot medium fill the same volume. This region corresponds to 
the disk-halo interface described above. In what follows, we 
shall consider the disk-halo interface as the height at which 
the two media have the same volume-filling factor. 

As the warm neutral and ionized media dominate 
different regions, it gives the appearance that each of these 
phases constitute different layers. However, both media are 
mixed and have different scale-heights. The region domi- 
nated by the ionized medium is known as the Reynold's 
layer (Reynolds, 1987). 

The filling factors of the different phases observed in this 
study are similar to those estimated by Spitzer (1990) and 



Ferriere (1995). These authors estimated a volume-filling 
factor of the hot gas in the Galactic disk of ~ 0.2. In addi- 
tion, Ferriere (1995) estimated volume-filling factors for the 
cold and neutral gases in the disk similar to those shown in 
Figure ^ and discussed above. 

Ferriere (1998) estimated the variation of the volume- 
filling factor of the hot gas with z. She finds that, at the 
solar circle, the volume occupation of the hot gas increases 
from 0.15 at z = pc to 0.23 at z — 200 pc and decreases 
gradually for z > 200 pc (see Figure 12 in Ferriere, 1998). 
The latter result seems to be in confiict with the calculations 
reported in this paper. Such apparent confiict is resolved if 
one considers that the calculations carried out by Ferriere 
refiect the distribution of hot gas inside cavities of supernova 
remnants and superbubbles, whose distribution decreases 
with 2, while the present paper also considers hot gas that 
bouyantly rose from interacting supershells. Ferriere's model 
does not include that gas at all, and so does not model the 
high- 2 region fully. 

4.2 Large Scale Outflows and Local Outbursts 

The simulations show that the thick gas disk is fed with 
the hot gas coming from the thin disk through two major 
processes: large scale outfiows and local outbursts (known 
as chimneys). 

Large scale outfiows are the result of the buoyant expan- 
sion of hot disk gas concentrated in large reservoirs located 
on either side of the thin HI disk (as shown in the density 
profiles). These reservoirs are formed by the hot inner parts 
of isolated supernovae randomly distributed in the stellar 
disk. 

Such hot gas has enough energy to be held gravita- 
tionally and expands buoyantly upwards, generating large 
scale outfiows (Figure The ascending fiow triggers the 
growth of Rayleigh- Taylor instabilities as it interacts with 
the cooler, denser medium in the thick disk. Therefore, the 
hot gas acquires a finger-like structure with a mushroom 
cap on the tip of the finger. During the expansion of the 
finger, further instabilities develop, until the major fraction 
of the ascending gas cools and merges with the surrounding 
medium. 

The gas in the finger cools from the outside towards the 
centre, thus giving it the appearance of being enveloped by 
a thin sheet of cooler gas (Figures |^ and ^ at x — 300, y — 
100 pc). A structure called the "anchor", with properties 
similar to those described above, has been observed in the 
southern Galactic hemisphere and reported by Normandeau 
& Basu (1998). 

Chimneys result from correlated supernovae, which gen- 
erate superbubbles that acquire an elongated shape due to 
the stratification of the ISM in the z— direction. Superbub- 
bles are able to expand and blow holes in the disk, injecting 
high speed gas, which forces its way out through relatively 
narrow channels with widths of some 100 - 150 pc. They 
provide a connection between the thin disk and the upper 
parts of the Reynolds layer. 

One can identify mini-chimneys and chimneys from the 
break up of shells and supershells located at the upper parts 
of the thin disk. Mini-chimneys result from isolated super- 
novae that have had their shells broken and released to 
hot gas into the surrounding medium (the broken shell acts 
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as the walls of a small tunnel produced by the blow out), 
whereas chimneys are large tunnels of hot gas crossing the 
thick HI disk. These structures are well coUimated, with a 
maximum width varying between 50 pc for mini-chimneys 
and 150 pc for chimneys, and vortical extensions varying 
between few tens of parsec (mini-chimneys) and 500 pc or 
more for the chimneys. 

This is corroborated by the HI position maps, taken at 
different locations of the disk along lines of sight parallel to 
the plane, showing the presence of broken shells ("worms") 
wiggling around in the ^-direction. Face-on HI maps of the 
disk taken at \z\ ~ 250 pc show the presence of low density 
regions with elliptical or circular forms and having a large 
range of sizes. The major fraction of these holes have been 
produced by large cavities of hot gas resulting from clustered 
supernovae, and therefore by chimneys. It is hard to identify 
the holes resulting from mini-chimneys, suggesting that their 
maximum elongation is some 150 pc in the z-direction. 

Surveys of the Galactic disk (Heiles, 1984, 1990) and 
of M31 (Brinks & Bajaja, 1986) and M33 (Deul & Hartog, 
1990) have revealed the presence of roughly elliptical holes 
in the HI maps, with sizes varying between 40 and 1000 pc. 
Some holes show a clear shell structure, although for most 
holes the shells are not resolved. 

The majority of the shells have deviations in the local 
velocity field with expansion velocities of 10-30 km s~^. The 
estimated kinematic ages of these holes varies from 2.6 to 
30 Myr and the energy requirements are 10*^ to 10'"''' erg. 
A correlation between the holes in the HI maps with single 
HII regions and OB associations in these galaxies has been 
estabUshed by Deul & Hartog (1990). 

Heiles (1984) identified the presence of holes in the HI 
emission maps of the Galactic disk. In addition, sheets of HI 
gas with a vertical structure have been identified in these 
maps by Heiles and Koo et al. (1992). These vertical struc- 
tures have been associated with broken shells and supersholls 
in the Galaxy, and have been proposed as being the walls 
of chimneys by Tomisaka & Ikcuchi (1986) and Norman & 
Ikeuchi (1989). HI maps of NGC 891 show a series of low 
density channels, similar to the ones shown in this paper. 

Most recently, Normandcau et al. (1996) identified a su- 
perbubble in the Perseus arm underlying an association of 
nine O-type stars. The superbubble has been created by stel- 
lar winds from these stars (Basu et al.., 1999) - no supernova 
has been identified within the W4 region - and heis an elon- 
gated appearance resulting from the vertical stratification of 
the disk gas. Although in HI emission maps the region looks 
like a chimney, such a view is altered in the Ha emission 
maps of the region, which clearly shows the presence of a 
superbubble. 

4.3 Numerical Resolution 

Numerical resolution seems to be an issue in the simulations 
reported in this paper and in Paper II. The major constraint 
imposed by numerical resolution is the prevention of the 
formation of very small scale structures resulting from ther- 
mal instabilities and condensations due to radiative cooling. 
Furthermore, the amount of cooling may be substantially 
increased if the resolution is great enough to allow for the 
development of turbulent shear layers in regions of outflow 
and infall. 



The numerical approach adopted to minimize such 
problems combined the usage of a high-resolution scheme 
(the PPM) to advance the solution in time with the dy- 
namical refinement of the computational grid whenever re- 
quired. The latter was set up in the region located at 
— 150 < z < 250 pc (hereafter referred to disk-grid), whereas 
in the rest of the computational domain no dynamical re- 
finement was used - the maximum resolution being 10 pc. 
The resolutions in the disk-grid were 5, 2.5 and 1.25 pc (this 
was the smallest grid size allowed during the calculations). 

The PPM is a third-order Godunov scheme, which does 
not rely on artificial viscosity to spread sharp variations of 
the flow variables. These discontinuities are represented by 
sharp profiles, free of spurious oscillations, of the flow vari- 
ables and are treated as solutions to the Riemann problem. 



5 FUTURE EFFORTS AND FINAL REMARKS 

The major drawback of the simulations is the absence of 
magnetic fields that may constrain the upward motions 

of the disk gas. Understanding the disk-halo interaction 
requires simulations that couple together the major con- 
situents of the ISM in the disk and halo, i.e, gas, magnetic 
fields, heating and cooling. 

The presence of magnetic fields may provide a confine- 
ment effect over the disk gas, providing the scale height of 
its distribution is much larger than that of the disk gas. For 
an average density distribution such as that observed in the 
simulations, the scale height is of the order of Hg ^ 170 pc, 
which is comparable to the scale height of the thin compo- 
nent of the Galactic magnetic field - Hb ~ 170 pc (Thomp- 
son & Nelson, 1980). Tomisaka (1998), using a magnetic 
field decreasing as the square root of the gas density, shows 
that no confinement occurs for the superbubblcs in the disk. 
However, the presence of a thicker component of the mag- 
netic field, with a scale height of 1.2 kpc (Han & Qiao, 1994), 
confines the expansion of superbubbles, located at \z\ < 300 
pc, to more than 20 Myr. 

Furthermore, as has been shown here, such study must 
be three-dimensional and independent of the grid size. That 
is, a 3D grid with a thickness of only a few hundred par- 
sec generates an overpressured region in the disk, and thus 
the simulations evolve in a two-dimensional fashion rather 
than three-dimensional. The time delay between supernovae 
and superbubbles is so large that when the next supernova 
explodes, the previous one has already expanded to few hun- 
dred parsec in size, and therefore crossed the boundaries of 
the grid. The younger supernova has a large probability of 
occurring inside the old cavity and, as such, it expands in 
an evacuated medium with velocities larger than that of a 
supernova in phase II. 

Such simulations are difficult to carry out due to 
present-day computational facilities. However, efforts are 
under way, by the author, to provide a new model where 
magnetic fields come into play in a medium where the disk- 
halo structure is regulated by supernovae. 

Phenomena such as high and intermediate velocity 
clouds are regarded as products of the Galactic fountain 
and their dynamics, structure and evolution are described 
in paper II (Avillez, 1999b). In addition, some of the items 
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described in this paper require full discussion which will be 
presented in future papers in the series. 
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